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Abstract. This paper is the first part of a work which consists in proving 
the stabilization to zero of a fluid-solid system, in dimension 2 and 3. The 
considered system couples a deformable solid and a viscous incompressible 
fluid which satisfies the incompressible Navier-Stokes equations. By deforming 
itself, the solid can interact with the environing fluid and then move itself. The 
control functions represents nothing else than the deformation of the solid in 
its own frame of reference. We there prove that the velocities of the linearized 
system are stabilizable to zero with an arbitrary exponential decay rate, by a 
boundary deformation velocity which can be chosen in the form of a feedback 
operator. We then show that this boundary feedback operator can be obtained 
from an internal deformation of the solid which satisfies the linearized physical 
constraints that a self-propelled solid has to satisfy. 



1. Introduction. In this two-part work we are interested in the way a solid im- 
mersed in a viscous incompressible fluid (in dimension 2 or 3) can deform itself and 
then interact with the environing fluid in order to stabilize exponentially to zero the 
velocity of the fluid and also its own velocities. The domain occupied by the solid at 
time t is denoted by S{t). We assume that S{t) C O, where O is a bounded smooth 
domain. The fluid surrounding the solid occupies the domain O \ S{t) = -7^(i)- 



O = T{t) U S{t} c K- or R''. 
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1.1. Presentation of the model. The movement of the soHd in the inertial frame 
of reference is described through the time by a Lagrangian mapping denoted by Xs, 
so we have 



where the vector h{t) describes the position of the center of mass and R(t) is 
the rotation associated with the angular velocity of the solid. In dimension 3 the 
angular velocity is a vector field whereas it is only a scalar function in dimension 
2. However, can be immersed in and this scalar function can be read on 
the third component of a 3D-vector. More generally in this work, since all the 
calculations made in dimension 3 make sense in dimension 2, we will consider only 
vector fields of M.^ and matrix fields of M'^^'^. For instance, u> and R are related to 
each other through the following Cauchy problem 



where in dimension 2 we have uj = and wi = W2 = 0. 

The couple (ft,(i),R(<)) describes the position of the solid and is unknown, whereas 
the mapping X*{-, t) can be imposed. This latter represents the deformation of the 
solid in its own frame of reference and is considered as the control function on which 
we can act physically. When this Lagrangian mapping X*{-,t) is invcrtible, we can 
link to it an Eulerian velocity w* through the following Cauchy problem 



The fiuid flow is described by its velocity u and its pressure p which are assumed 
to satisfy the incompressible Navier-Stokes equations. For X* satisfying a set of 
hypotheses given further, the system which governs the dynamics between the fluid 
and the solid is the following 



S{t) ^ Xs{S{Q),t), t>0. 
The mapping can be decomposed as follows 

Xs{y,t) = h{t) + Ii{t)X*{y,t), yG5(0), 





X*{S{Q),t). 



— i^Au + {u ■ V)u + Vp = 0, 

ot 

div u — 0, 



X e T{t), 



t 6 (0,oo), 
t G (0,oo), 



(1) 
(2) 



M = 0, 

h'{t) + Lu{t) A (x - h{t)) + w{x, t), 



xedO, te(0,oo), 
xedS{t), t(E{0,oo) 



(3) 
(4) 




i G (0,oo), 



t G (0, oo). 



(5) 



(6) 



u(y,0) = wo(y), y G J'(O), /i'(0) = /ii G M^ w(0) = wo G ]R^ 



(7) 
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where 

S{t) = h{t)+Rit)X*{SiO),t), J^{t) = 0\S{t), (8) 

and where the velocity w is defined by the foUowing ehange of frame 

w{x,t) = R{t) w* {R{tf{x - h{t)),t) , xeS{t). (9) 

Without loss of generality, we assume that h{0) = 0, for a sake of simplicity. The 
symbol A denotes the cross product in R^. The linear map w A ■ can be represented 
by the matrix S(a;). In equations (5) and (6), the mass of the solid M is constant, 
whereas the inertia moment depends a priori on time. In dimension 2 the inertia 
moment is a scalar function which can be read on the inertia matrix given by 

I{t) = (^J^^^psix,t)\x-h{t)\''dx^l^.. 

In dimension 3 it is a tensor written as 

I{t) = / psix,t){\x-h{t)\HM3 -{x-h{t))(g){x-h{t)))dx. 
Js{t) 

The quantity ps denotes the density of the solid, and obeys the principle of mass 
conservation 

= det(VX;L)) ' '''^'^^ 

where VX5 denotes the Jacobian matrix of the mapping ■ For a sake of simplicity 
we assume that the solid is homogeneous at time t = 0: 

psiy,0) = ps>0. 

In system (l)-(9), ly is the kinematic viscosity of the fluid and the normalized 
vector n is the normal at dS(t) exterior to J-{t). It is a coupled system between the 
incompressible Navier-Stokes equations and the differential equations (5)-(6) given 
by the Newton's laws. The coupling is in particular made in the fluid-structure 
interface, through the equality of velocities (4) and through the Cauchy stress tensor 
given by 



cr(w,p) = 2vD{u) -pld^v {Vu + [Vuf^ - p 



Id. 



Indeed, the Dirichlet condition (4) partially imposed by the deformation of the solid 
(through the velocity w) influences the behavior of the fluid whose the response is 
the quantity (j{u,p)n in the fluid-solid interface. It represents the force that the 
fluid applies on the solid, and then it determines the global dynamics of the solid 
(through equations (5) et (6)) and thus its position. 

The problem is the following: What is the deformation X* of the solid we have to 
impose in order to stabilize the environing fluid and thus induce a behavior of the 
fluid which stabilizes the velocities of the solid? We shall assume a set of hypotheses 
on the control function X* , that we state as follows: 

HI: For aU t > 0, X*{-,t) is a C^-diffeomorphism from S{0) onto S*{t). 
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H2: In order to respect the incompressibility condition given by (2), the volume 
of the whole solid has to be preserved through the time. That is equivalent 
to say that 

BX* 

— — • (cofVX*) ndr = 0. (10) 

05(0) Ot 

H3: The linear momentum of the solid is preserved through the time, which 
leads to 



P5(y,0)X*(y,0dy = 0. (11) 

/5(0) 

H4: The angular momentum of the solid is preserved through the time, which 
leads to 

BX* 

ps{y,0)X*{y,t)^—-{y,t)Ay - 0. (12) 

5(0) Ot 

Imposing constraints (11) and (12) enables us to get the two following constraints 
on the undulatory velocity w 

ps{x,t)w{x,t)dy = 0, (13) 

5(t) 

psix,t){x ~ h{t)) Awix,t)dy = 0. (14) 

5(t) 

As equations (5) and (6) are written, the constraints (13) and (14) are implicitly 
satisfied in system (l)-(9). Hypotheses H3 and H4 are made to guarantee the 
self-propelled nature of the motion of the solid, that means no other help than its 
own deformation enables it to interact and to move in the surrounding fluid. 
The existence of global-in-time strong solutions for system (l)-(9) has been studied 
in [19] in dimension 2 and more recently in [5] in dimension 3. In particular, this 
existence in dimension 3 is conditioned by the smallness of the data, that is to say 
the initial condition (uq, /ii, wq) and the displacement of the solid X* — Ids (in some 
Sobolev spaces). 

1.2. The linearized problem. For the full nonlinear system (l)-(9), the equa- 
tions are written in the Eulerian configuration, and thus we are lead to think that 
the Eulerian velocity w* is the more suitable quantity to be chosen as a control 
function (instead of X*). But such a mapping is defined on the domain S*{t), 
which is itself defined by X*{-,t). Moreover, the study of such a nonlinear system 
is based on the preliminary study of the corresponding linearized system which is 

BU 



Bt 



div cr(?7, P) = 0, in J"(0) X (0,oo), (15) 

div U ^Q, in J"(0) X (0, oo), (16) 



U^Q, in 90 X (0,oo), (17) 

U^H\t) + n{t)Ay + ay.t), 2/6 95(0), te(0,c»), (18) 

MH"{t) = - I cr(C/,P)ndr, <e(0,oo), (19) 

JdS 

Ion'{t) = - I yAa{U,P)ndr, ie(0,oo), (20) 

JdS 
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U{y, 0) = m{y), y G ^(0), i/'(0) = /ii e M^, 1^(0) = ujo e R^, (21) 
and where the more suitable control to be chosen is the function C, related to the 

, ■ dx* 

Lagrangian velocity — - — by 
at 



C = e 



At 



dX* 



dt \ds ' 

Notice that the constraints (10) and (12) are nonlinear with respect to the mapping 
X* . We linearize them when we consider the linear system (15)-(21). For this 
linear system, the constraint induced by Hypothesis HI can be relaxed, since we 
only consider mappings X*{-,t) continuous in time and such that X*{-,Q) = Ids- 
Thus the notion of admissible control for this linearized is made precise in Definition 
2.1. 

1.3. The main result and the strategy. The main result of this first part is 
Theorem 6.1, which implies the following one: 

Theorem 1.1. Assume that dist(3C'; 5(0)) > 0. For all (uo,ft.i,a;o) satisfying 
Uq G H^(J^(0)) and the following compatibility conditions 

div uo = in J-{0), 

uo = on do, 

Uq ^ hi + u)q a y on dS{0), 

system (15)~(21) is stabilizable with an arbitrary exponential decay rate A > 0, that 
is to say that for aU A > there exists a boundary control ( e L^(0, oo; H^/^(95)) fl 
H-'^(0, oo; H-^^^idS)) and a positive constant depending only on {uq, hi,ujo) such 
that for all t > the solution (U, H' , ^l) of system (15)-(21) satisfies 

||([/(-,t),i7'(t),O(t)||L2(^(0))xR3^R3 < Ccxp(-At). 

For proving this theorem we study the system that has to satisfy the functions 

Ui-, t) = e^*[/(-, t), P{-, t) ^ e^*F(-, t), H'{t) ^ e^*H{t), n{t) e^'n{t), 

and the goal is then to prove that there exists a control C such that this system 
admits a solution (U,H',Cl) bounded in some infinite time horizon space. The 
strategy we follow is globally the same as the one used in [17], at least for the 
linearized problem. It first consists in rewriting the full nonlinear system in space 
domains which do not depend on time anymore, by using a change of variables and 
a change of unknowns. Then we can make appear all the nonlinearities (specially 
those which are due to the variations of the geometry through the time) and we 
can set properly the linearized system. The second step of the proof consists in 
formulating the linearized system in terms of operators where the pressure is actu- 
ally eliminated and encodes a mass-added effect. This writing enables us to define 
an analytic semigroup of contraction generated by an operator which presents in- 
teresting spectral properties: Indeed the unstable modes that we have to stabilize 
are actually countable and in finite number. Besides, we prove by a unique contin- 
uation argument the approximate controllability of this linearized system. Thus, 
in order to define the boundary control that stabilizes the full linearized system, it 
is sufficient to consider a finite-dimension linear system for which the approximate 



^ In the following the symbol C will denotes a generic positive constant which does not depend 
on time or on the unknowns. 
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controllability is equivalent to the feedback stabilizability. The aforementioned con- 
trol can be defined on a finite-dimension space in a feedback operator form, what 
will be useful for proving the stabilization of the full nonlinear system. 
With regards to the methods, a novelty is the means provided in a last section 
which enables us to define from a boundary control an internal deformation satis- 
fying the linearized constraints. This result too will be useful for the definition of a 
deformation of the solid - which has to satisfy the nonlinear constraints - that stabi- 
lizes the nonlinear system: Considering a deformation which satisfies in a first time 
the linearized constraints is not necessarily from a mathematical point of view, but 
the method with which we obtain it will be important for defining a deformation 
satisfying the nonlinear constraints; Besides, from a physical point of view, consid- 
ering for the linearized system a deformation which satisfies the linear constraints 
is relevant since it ensures the conservation of the momenta for the whole fluid-solid 
linear system. 

The idea of considering first the linearized problem relies on the fact that for small 
perturbations (that is to say for small initial conditions uq, hi and wq) the behavior 
of the nonlinear system is close to the one of the linearized system. 

Thus the same statement will be proven in the second part of this work for the 
unknowns of the full nonlinear system (l)-(9). The result is nonintuitive: It says 
somewhat that all the fluid in which the solid swims can be stabilized just by the 
help of this swimmer, at an intermediate Reynolds number. This kind of problem 
has been investigated in [10, 11] for instance, where it is considered other types of 
fluid-swimmer systems. The same kind of purpose has been also investigated at 
a low Reynolds number in [18] and more recently in [12], for the Stokes system, 
and in [3, 4] in the case of a perfect fluid. The control of the motion of a boat at 
a high Reynolds number has been recently studied in [8]. Besides, the same kind 
of techniques that we use in this work have been used for other coupled systems 
involving the incompressible Navier-Stokes equations; Let us cite [1] and [17] for 
instance, where the stabilization of other fluid-structure problems is proven. 

1.4. Plan. Deflnitions and notation are given in section 2. Next in section 3 we 
explain how to rewrite the main system in space domains which do not depend 
on time. Then we will be in position to linearize the system so obtained. The 
linearized system is studied in section 4 where it is rewritten through an operator 
formulation for which we prove useful properties. The approximate controllability 
of this linearized system is proven in section 5. It leads to the main result of this 
paper, that is to say the feedback stabilization of the linearized system in section 
6. Finally in section 7 we provide a means to recover an internal deformation 
of the solid from a feedback boundary control which lives only in the fluid-solid 
interface. This final section is the transition to the second part of this work where 
the considered control is an internal deformation of the solid. 

2. Definitions and notation. We denote by = T{0) the domain occupied by 
the fluid at time t = 0, and by 5 = 5(0) the domain occupied by the solid at t = 0. 
We assume that S is smooth enough. We set 

Let us introduce some functional spaces. We use the notation H'*(i7) = [II*(ri)]'' or 
[II'*(il)]'* , for some integer k, for all bounded domain ft of or K.^. We classically 
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define 

V°(J-) = {0eL2(J-) I div = in J", -71 = 01130}, 

Vi(J-) = {0gH1(J') |div<?i = Oin 0-n = OonaO}, 

U'-HQl) = L2(0,^;H2(^))nHi(G,cx);L2(^)), 

Let us keep in mind the continuous embedding: 

H^'^Ql) ^ L°°(0,oo;Hi(^)). 

We finally set the spaces dealing with compatibility conditions 

= {(mo,^i,wo) e V°(J") X R3 X I uo = /ii +u;o Ay on 55} , 

Hj^ = {{uo,hi,ujo) eVl.iJ') xR^ xR^ \ uo^ hi+ujo Ay ondS} . 

For more simplicity, wc assume that the density ps at time t = is constant with 
respect to the space: 

Ps{y,0) = P5 > 0. 

We assume without loss of generality that /iq = 0. This implies in particular 



ydy = 0. 

Let us state the conditions we shall assume on mappings chosen as control functions. 

Definition 2.1. A deformation X* is said admissible for the linear system (15)-(21) 
if 



and if for alH > it satisfies the following hypotheses 

dX* , , 



dS 



dt 



f X*{y,t)dy = 
Js 



oo;Hi(5)) 




0, 


(22) 


0, 


(23) 


0. 


(24) 



dX* , , , 
y A -^(y,t)dy 

3. Rewriting in cylindrical domains and linearization. System (l)-(9) is 
strongly coupled, in particular from a geometrical point of view. A first work then 
consists in rewriting it in space domains which do not depend on time anymore. 
For that, we extend to the fluid domain the change of variables Xs induced by the 
deformation of the solid. 

3.1. The change of variables. The aim is to define a mapping 

X{-,t) : T — >T{t) 

which has to be for alH > a C^'^-diffcomorphism from T onto J-'(<), and which has 
for alH > the following properties 

detVX(?/,t) = 1, y G J", 
X{y,t) = y, yedO, 
X{y,t)=Xs{y,t), yedS. 

For that, let us first define an intermediate extension. We extend the mapping X* 
to F as follows: 
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Proposition 1. Let us assume that the mapping X* satisfies the hypothesis H2, 
that is to say for all t > we have the condition 



f dX* 

/ — — • (cofVX*)ndr = 0. 
Jds (^t 




dX* 

Let us also assume that the function [y, t) ^ e^* is small enough in L'^(0, oo; H'^(iS))n 

Hi(0,oo;Hi(5)). Then there exists a mapping X* e L°°(0, cx); H3(J-))nWi'°°(0, cx); H^J^)) 
satisfying 

in T (0, oo), 
on dS X (0, oo), 
on do X (0, oo). 

Proof. We give the proof of this result in Part II of this work; Indeed, the technical 
details of the proof and the estimates which follow are not useful in this Part I, but 
will be essential for proving the main result of Part II. In this paper we just need 
to know that the extension provided by this proposition exists. □ 

The mapping Xg^-., t) is obtained from i) by composing it to the left by the 
rigid displacement 

X^\-,t) : X* ^ hit) + 'R{t)x*. 

We cannot do the same thing for obtaining the mapping X from X because of the 
boundary condition on dO that has to be preserved. Thus we define an extension 
of X^{-,t) to the whole domain. For that we can use the same process which has 



been introduced in [20], and thus construct X (■,i) which satisfies for alH > the 
following properties: 

deWx"'i-,t) = 1 inR^orM^, 
X''(.,t)=X^(.,t) ond{X*{S,t)), 
x"{-,t)^ldao on do. 

Here again in a sake of simplicity we do not give more detail on this construction, 
since it is useless in this Part I, and it is completely detailed in Part II. So we define 
X as follows 

X{;t) = x'''{;t)oX*{.,t), 

and we denote by Y{-, t) its inverse for all t > 0. 

Remark 1. The definition of such a change of variables is also done in [19] or in 
[5]. The change of variables we have constructed here has the same properties as 
the ones constructed in these articles. But the way we proceed here is not the same: 
In [19] the mapping X is constructed by extending the Eulerian velocity w to the 
fluid part, but this means is not suitable in a framework where the regularity of 
the deformation of the solid is limited; Indeed the velocity w{-, t) is defined through 
i) which is itself defined on the domain X*{S, t). Concerning the means used in 
[5], the mapping X is directly defined without extending X* first, but the existence 
of X is obtained only locally in time. Here the hypothesis we make is the smallness 
of e^*^^ in an infinite time horizon space. 
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In the following, we will denote 

X{y, t) = R(t)^(X(y, t) - h{t)), (y, t) e T x (0, 

Y{x, t) = Y{h{t) + R(t)x, t)n {x, t) e [j F{t) x {t). 

t>o 

3.2. Rewriting of the main system in a cylindrical domain. Let us transform 
system (1)~(9) into a system which deals with non-depending time domains. For 
that we use the change of variables defined above in order to set the following change 
of unknowns 

u{y, t) = 'R.{t)^u{X{y, t),t), u{x, t) = R{t)u{Y{x, t), t), 

p{y, t) = p{X{y, t),t), p{x, t) = p{Y{x, t), t), 



for X G ^{t) and y £ J^, and 

h'{t) ^ Il{tfh'(t), Cj{t) = Ti(tfuj{t). 

This is the same change of unknowns which is used in [5] for rewriting the main 
system in cylindrical domains. For a sake of clarity, let us also define the mappings 

X{y, t) = R(t)^(X(y, t) - h{t)), Y{i, t) = Y{h{t) + R(t).i, t). 

Then, like in [5], system (l)-(9) whose the unknowns are h' is rewritten in 
the cylindrical domain J- x (0, cxd) as the following system whose the unknowns are 
{u,p,h' ,u)): 

du ~ 

lyLii + M(u,h',uj) + Nu + uj(t) A it + Gp = 0, in J" x (0, oo), 

ot 

div u = g, in 7^ x (0, cxo), 

■u = 0, in do X (0, oo), 

dX* 

u = h'{t) + Co{t) A X*{y,t) + ^{y,t), (y, t) £ 95 x (0, oo), 

Mh"{t) = - / a{u,p)'^Y{XfndT - MCo{t) A h'{t), t £ (0, oo) 

JdS 

I*{t)C'{t) = - [ X*{y,t) A (~a{u,p)VY{Xfn) AT 
JdS ^ ' 

-~r'{t)Cj{t) + I*{t)Cj{t) A (L{t), t e (0, oo) 

u{y,Q)^uo{y). y£T, h'{Q)^hi£R\ lD(0) = e M^, 
where [-Ji specifies the i-th component of a vector 
[l.u],{y,t) = [Vu{y,t)l^Y{X{v,t),t)], + \I^u,{v,t) : (vfvf^) (l(y,i),0, 

M(u, h', Co){y, t) = -Vu{y, t)yY{X{y, t),t) (h'{t) + ib{t) A X{y, t) + ^{y, t) j 

Nu{y,t) = Vu{y,t)VY{X{y,t),t)u{y,t), 
Gp{y,t) ^\/Y{X{y,t),tf\/p{y,t), 

d{u,p){y, t)^v (Vu(y, t)VY{X{y, t),t) + Vf(l(y, t), tfVu{y, t)^) - p{y, t)h 
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and 

g{y,t) = trace ( Vw(y, i) ( Ir3 -Vy(X(y,t),t 



Notice that we have 

cof(VX) dct{\7X)VY{Xf = VY{Xf, 

and so from the Piola identity we can actually express this nonhomogeneous diver- 
gence term as g ~ div G, where 

Giy,t) = (lM3-Vy(l(y,i),i))"(y,i)- 
For A > 0, we make the following change of unknowns: 

r) X* 

u = e^*u, p = e^'p, h' = e^*h', w = e^*w, C = e^*-^. (25) 

The idea of this second change of unknowns is the following: If we find a control X* 
such that the quadruplet {u,p, h' , uj) is bounded in some infinite-time horizon space, 
then the intermediate unknowns {u^p, h',uj) will be stabilized with an exponential 
decay rate, and it will be sufficient to deduce from that the same property for the 
actual unknowns {u,p, h' , uj), in a functional framework which will be made precise 
in the second part of this work. 

The system satisfied by {u,p, h',uj) is then transformed into 

^-Au-i/Lu + e-^*M(?i,/i',w) + e-^*NM + Gp = F, y e T, te(0,oo), 

ot 

div u = div G, y G J-', t £ (0, oo), 

M = 0, y&dO, te(0,oo), 

u = h'{t)+CJ{t)^X*{y,t) + (^{y,t), yedS, te(0,oo), 

Mh" [t) - \Mh' = - I a{u,p)VY{XY'ndV + Fm, te{Q,oo) 
Jas 

r {t)Cj' {t) - \I* {t)uj{t) ^ - f X* AU{u,p)VY{Xfn)dV + Fi, t € (0, oo) 

Jds ^ ' 

u{y,Q) = uo{y), y gT, /i'(0) = /ii e M^, u;(0)=woeM3, 

with 

F = -e~^*(2; A u, Fm = -Me'^^Cb A h' , Fj = -F'uj + e^^^ Tuj A uj. 

In order to consider the linearized problem associated with this nonlinear system, 
we rewrite this latter as follows 
Ou 

— \u — lyAii + Vp = F{u,p, h' ,Cj), in J" x (0, oo), 

at 

div u = div G(u), in x (0, oo), 

w = 0, on do X (0, oo). 



u 



= h'{t) + uj{t) Ay + C, + W{uj), {y, t) e dS x (0, oo). 
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Mh"-\Mh' = - (j{u,p)ndr + FM{u,p,h',uj), in(0,oo), 
Jas 



IqCj' [t) ~ XIqOJ — — / y A (T({t,p)ndr + w), in(0, oo), 

JdS 

u{y,Q) = uo{y), ijeF, h'{Q)^h^^WL\ Cj{Q) ^ lu^ eB? , 

with 

F{u,p,h',Cj) = i/(L-A)M-e-^*M(w,A',w)-e-^*Nit-(G-V)p-e-^*wAu, 
G{u) = (lR3-Vy(l(y,i),i))^i, 
W{lo) = ljA{X*- Id) , 
FM{u,p,h\u}) = -Me^^^Lo Ah'{t) 

-V j (vu ivY{X) - Irs) + (VY{X) - Irs)^ Vu^^ Vf (X)^ndr 

- / <j{u,p){vY{X)^Ik.Y ndV, 
JdS ^ ' 

Fi{u,p,u}) = -{I* - h)uj' + \{r - h)uj - r'Co + e-^^rCo ALU 



-V j yA [vii [VYix) - Ir3 j + (vy(x) - Ir3)^ Vu^ j vr(x)^ ndr 



y A (^a{u,p){VY{X) -l^.fn^ dV 
+ I (X* - Id) A (w*{u,p)VY{X fn^ d 



3.3. Linearization. 

3.3.1. Statement. We linearize system (26)-(26) around {u,p, h\ w, X*) = (0, 0, 0, 0, Id), 
and wc obtain formally 

dU ' . . 

— ~XU-iyAU + VP = 0, in J" X (0,cx)), 
ot 

div U = 0, in X (0, cxj), 

C/ = 0, in90x(O,oo), 
U^H'it) + n{t)Ay + C{y,t), y e dS, t e {0,oo), 

MH"{t) - \MH'{t) = - / a{lJ, P)ndr, t e (0, oo), 

lon'it) - XloClit) = - [ yAa{U,P)ndr, te(0,oo), 

?7(2;,0) = uo(2/), yG^, i?'(0) = /ii G M^ f2(0) = G K^, 

with 

a(;7, P)(zj, t) = 2,yD{U){y, t) - P{y, t)ld. 

Note that in this linear system the control (initially chosen as the deformation of 

dX* 

the solid) appears only on the boundary dS through the function ^ = e'*'* 



dt 
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Thus the problem of stabilization of this linear system is reduced to a problem of 
boundary stabilization. 

3.3.2. Constraints adapted to the linear system. Let us observe that the constraints 
(10) and (12) are nonlinear, with respect to the mapping X* . These constraints are 
written 

r 3x* r 3x* 

Jos ~dr ■ *)^) = 0, j^X*A -^dy = 0. 

The linearization of this equalities leads us to consider the following constraints 

dX* f dX* 

— - • ndr = 0, y^ —-dy = 0. 

ds dt Js dt 

The constraint (11) is already linear. This leads us to consider the constraints 
stated in Definition 2.1 for controlling the linear system (15)-(21). 

4. Operator formulation: Definition of an analytic semigroup. In this sec- 
tion we study the following linear system, which is nothing else than system (15)- 
(21) with a resolvent term: 

du 

— Am — div (t(u, v) = 0, in , 

dt ^ ' 

div u — 0, in , 

u = 0, in 9( 

u = h'{t)+uj{t)Ay + C{y,t), ye, 



(0,oo), 


(26) 


(0,oo), 


(27) 


(0,^), 


(28) 


t e (o,oo). 


(29) 


in (0, oo). 


(30) 


in (0, oo). 


(31) 


a;(0) = Wo e R^. 


(32) 



Mh" - \Mh' = - / a{u,p)ndT, 
JdS 

Iquj' — XlctOJ = y A a{u,p)ndr, 

JdS ' 

u{y,0) = uo{y),yeJ^, h'{0) = hieM.\ 

4.1. Introduction of some operators. For what follows, we need to introduce 
some operators. Let us first recall the notation of the stress tensor of some vector 
field u 

D{u) = i (Vu + Vu^) . 
and let us denote the Hessian matrix operator as 

H = V^. 

For h' G R"^ and £ R'^, we define N{h') and N{uj) as being the respective solutions 
q and q of the Neumann problems 

Ag = in ^^h' -n on dS, -^^ = on dO, 
on on 

dq dq 
Aq = in J". — = (uj Ay) ■ n on dS, — — = on dO. 
on On 

For ip £ Ti^/^{dS), we define Laip = w as being the solution of the Stokes problem 

— z^Aw + VtA = in J^, div w = in J-", 

w = on do, w = If on dS. 
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Similarly we define L^Lp - 
— i^Aw 



w as being the solution of 

- VtA = in J^, div w = in J-", 

w = on do. 



w ~ if Ay on dS. 



The operators Lq and Lq are called lifting operators. We also define the following 
integration operators 



Cip = / ipndT, 
'as 



Cif = / y A i^Tidr. 

Ids 



We denote by P : L^(J^) t-s- V^j(J^) the so-called Leray or Helmholtz operator, which 
is the orthogonal projection induced by the decomposition 

Then we denote in V° the classical Stokes operator by 

Ao = z/PA, 

with domain D{Ao) = H2(J-) n Hi(J-) n (J"). 

4.2. Operator formulation. Let us first consider system (26)-(32) only when 
C = 0; In that case we denote by {v,p, h' , uj) the unknowns. But this system can be 
transformed into a system whose unknowns are only {v, h' ,lu). Indeed, by following 
the method which is used in [16] or [17] for instance, the pressure p can be eliminated 
in the equations (26), (30) and (31). By this means we obtain that p can be written 



p = t: 



at' 



where tt is solution of the following Neumann problem 

dnit) 



A7r(i) = in J", 



9?- 



= vFAvit) ■ n on dT, 



and q is solution of this other Neumann problem which involves the boundary 
conditions 



Aq{t) =^ in T, 



^''^'^ -{h'{t)+uit)Ay)-nondS, ^""^'^ 



on do. 



dn dn 
Moreover, Vp can be expressed through a lifting; More precisely, we have 

Vp = {-Ao)PLa{h' +UJ Ay) in T. 

Thus, we can split system (26)-(32) into two systems, one satisfied by {Vv,h',u!), 
and the other one satisfied by (Id — ¥)v. Explicitly, by denoting V = {Vv, h',ijj)'^, 
system (26)-(32) can be rewritten (for C = 0) as follows 

iMa+Madd)V' = AV + XMoV, (33) 

(Id-P)^; = {ld-P)(^Lo{h') + Lo{uj)) , (34) 



with 



Id 
Mid 
/o 



^add 



■ 

CN CN 
CN CN 
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and 

A = C{-2iyD + NAo) -CHN -CHN 
C{-2vD + NAo) -CHN -CHN 

The operator M.add represents a mass-added effect. We are going to see that it 
contributes to making the "effective mass operator" - that is to say Mq + M^add - 
self-adjoint and positive. 

4.3. Main properties of the operator A. Let us set 
Lemma 4.1. M is self-adjoint and positive. 

Proof. Observe that Mq is self-adjoint and positive. Then it is sufficient to show 
that Marffi is self-adjoint and non-negative. Let us begin with noticing that CN is 
self-adjoint. Indeed, if qi and q2 denote respectively N{h[) and N{h2), by using 
twice the Green formula wc get 

/^Vgi-Vgzdy = / {h[ ■ n)q2dT = [ {h'^ ■ n)qidr 
Jas JdS 
h'^-CNih'^) = h'2-CN{h[). 

We can also see that CN is self-adjoint. If qi and (j2 denote respectively iV(cJi) and 
N{uj2), by using twice the Green formula we get 

Jjr Vgi • Vg2dy = / oji ■ {y A n)q2dT = / uj2 ■ {y A n)qidT 
JdS „ . JdS 

= UJi-CN(lU2) — U)2-CN{uJi). 

Likewise, let us show that (CN)'^ ~ CN. First, we denote q ~ N{uj), and by using 
twice the Green formula, we have 



CN{uj) = I Vgdy = / ^ydT 



T , JdS 



y^y{nAu})dT = / y ® {y A n)dr ] uj. 
dS \JdS / 

On the other hand, if we denote q = CN{h'), let us notice that q = h' -y. And then 

CN{h') = [ {yh'){yAn)dT = (f {y A n) ® ydT] h' . 

JdS \Jas J 

Then we can conclude that M.add is self-adjoint. In order to prove that M.add is 

non-negative, let us see that the corresponding quadratic term can be written 

V^MaddV = J^\yq + Vg|' dy > 

for V = {Pv,h',u})^, q = N{h') and q N{uj). □ 
In the following, we will denote 

A = (Mo+M,dd)"'A. 
Proposition 2. The domain of A is 

D{A) = n (H2(J-) X X ^ 

and A = is self-adjoint. 
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Proof. Let Vi = (wi, /li, wi)^ and V2 = (f2, /I'a, ^2)^ lie in D{A). We set {Fi,Fm, , Fi, Y 
(AM - K)Vi and {F2,Fm2,Fi^)'^ = (AM - 4)^2, that is to say that for i e {1,2} 
we have 

Xvi - div a {vi, Pi) = Fi, in J", 

div Vi = 0, in J^, 

= 0, in do, 

Vi = h[{t)+uj,(t) Ay, ye 35, 

AM/i^(t) = - / a(«,,p,)ndr + ^^m,, 

XIoUJi{t) = - y A<7{vi,pi)ndT + Fj^, 
JdS 



with 



Pi = N{Anv,) - \N{h'^ +oJ^^y)■ 

We calculate 

(V2; (-Fl,FA/i,i^7j'^)L2(jr)xR3xR3 = X{Vl , V2) 1,2 (jr) ~ j Vl ■ dw Cr(ui,Pl) 

+/12 ■ / CT(-(;i,pi)ndr + CJ2 • / y A CT(wi,pi)ndr, 

JdS JdS 

and by integration by parts we get 

Vi • div cr(wi,pi) + /12 ■ / cr{vi,pi)ndT + LU2 ■ / y A (T(ui,pi)ndr 

J^' JdS JdS 

V2 ■ a{vi,pi)ndT + 2iy D{vi) : D{v2) 

dS Js 

+/12 • / <7{vi,pi)ndT + UJ2 ■ j y A cr(ui,pi)ndr 

JdS JdS 



= 2v / D{vi) : D{v2). 
Js 

Then, by swapping the roles of Vi and V2 , it is easy to see that 

(V2; (i^i, i^A/i, i^/J^)L2(jr)xR3xR3 = {iF2, Fhh, ^h)'^ ^i)l2 (J^) xR3 xR3 . 

It shows that AM — A is self-adjoint. Since M is self-adjoint, the proof is complete. 

□ 

Proposition 3. The resolvent of A is compact. 

Proof For F = {F, Fm , Fif e L2(0, 00; L2( J")) x L2(0,oo;M3) x L2(0,oo;R3), we 
consider the system 

M (Aid -A)V = F, 
where V = (w, /i',a;)^ is the unknown. This system can be rewritten as 
Xv — div a{v,p) ~ F, in J-, 

div v = 0, in 

v^O, in do, 

v = h'{t) +iu{t) Ay, yedS, 
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XMh'{t) = - a{v,p)ndT + Fm, 

JdS 



XIoLo{t) 



y A a{v,p)ndT + F/ 



ds 



with 



p ^ N{Aov) - XN{h' + UJ Ay). 

By the same kind of ealeulations as the ones made in the proof of Proposition 2, 
this problem is equivalent to the following variational problem: 

Find V eV^xR^ X such that a{V, W) = 1{W) (35) 

with V = {v, h', oj)^, W ^ (w, k', a)^, and 



a(y, W) 



liW) = 



A / V-W + Mh' ■ k' + louj ■ a] + 2iy D{v) : D{w) 



F -w^Fu-k' ^Fi-a. 



By choosing A = 1, we use the Lax-Milgram theorem and prove that problem (35) 
has a unique solution. Thus M (Aid — A) is invertible, and since M is positive, 
Aid — .4 is also invertible. □ 

The results of the two last propositions yield the following theorem. 

Theorem 4.2. The operator {A, D{A)) is the infinitesimal generator of an analytic 
semigroup on V^{J-') x R"^ x Rp , and the resolvent of A is compact. 

4.4. Abstract formulation of the control problem. 

Proposition 4. The triplet {u,h',Lu) is solution of system (26)-(32) if and only if 
U = {¥u,h',uj)'^ and (Id — P)u satisfy the operator formulation 

U' = AxU + BxC, (36) 



(Id - P) M 

with Ax=A + AM-^Mo, Bx = 



(Id-P)(Lo(/i') + Lo(t^)) . (37) 

"^Ba = Ba; and 

(Aid - ^o) Lo 



Proof. Let us first recall that C must obey an incompressibility constraint given by 

C • ndr = 0. 

dS 

Thus we can formally extend C in the whole domain O while assuming that div C = 
in J-. That is why we can make the control appear only in the equation (36), the 
one which deals with Pu. Let us recall the linear system (26)-(32): 



— Xu — div a(u, v) 

dt ^ ' 



0, \n F X (0, oo), 
0, in J" X (0, oo). 



u = 0, in do X (0, oo), 

u^h'{t)+Lo{t)Ay + C{y,t), yedS, te(0,oo). 



STABILIZATION OF A FLUID-SOLID SYSTEM: PART I 



17 



Mh" - \Mh' 



as 



a{u,p)ndT, in(0,oo), 



IqUj' — \IqLO ~ — f y A (T{u,p)ndr, in(0,oo), 



OS 

The boundary condition on dS can be tackled by using a lifting method, like in 
[5] for instance. It consists in splitting the velocity u = v + w and the pressure 
p = g + TT. so that we have 

-i^Aw + Vtt 0, in J", 
div w = 0, in 

w = C, on dS, 
w = 0, on do, 



and 



dv dw 

— Xv — div a(v, q) = — - — h Aw, in x (0, oo), 

at ot 

div u = 0, in X (0, oo). 



w = 0, in aO X (0,oo), 

V = h'{t) +uj{t) Ay, y^dS, iG(0,oo), 



Mh"{t) - \Mh\t) = - / (j{v,q)nAT - / cr(w,7r)ndr, in (0,oo), 

JdS JdS 

IqU)' {t) — \I^uj{t) = — I yAa{v,q)ndT— / y A fT(w, 7r)ndr, in<£(0, oo), 
Jas JdS 

v{y,0) - uo{y)~w{y,0),yeJ^, h'{0)=hieM.\ w(0) = c^o e R^. 
This system can be formulated as follows 

MV^' = AF + AMoF + BiC + Ba,oC, 
{ld-F)u = (Id-P) (io(/i')+^H) , 

with V = (Pw, /i', w), C = II and 



-La 





ALo 

C{-2vD + NAo)Lo 
C{-2vD + A^^oi^o 



Notice that M"^Ba,i = Ba.i- The Duhamel's formula gives 

V(t) = e*'^MC^(0)+BiC(-,0))+^ e'*-")-^^ (M-iBA,oC + BiC)ds, 
and an integration by parts leads to 



iCds 



8iC(-,t)-e*-^^BiC(-,0)+ / e^'-'^^-^MABiCds 
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Notice that MqBi = and that 



Thus wc have 



Uit) 
U'it) 



C/(0) 



-C{-2vD + NAo)Lo 
-C{-2vD + N Ao)Lq 



g(t-s)^.j^-i (jg^^^ ^ ABi) Cds, 



AxU{t) + BxC,{-.t). 



and finally system (26)-(32) can be expressed formally in the form given by (36)- 
(37). □ 

4.5. Regularity of solutions of the linearized system. 

Proposition 5. For T > 0, {uQ,hi,UQ) G Hj^ and ( & L'^{0,oo;U^/'^{dS)) Ci 
H^(0, cxd; H^/^(9<5)), system (15)-(21) admits a unique solution {U,H\fl) such that 

Ueli^^^Q^), e h1(0,T;R3), n<Ell\0,T;m.^). 

Proof. With regards to the formulation (36)-(37) and the properties of the operator 
A, this result can be deduced from Proposition 3.3 of [21]. □ 

5. Approximate controllability of the homogeneous linear system. In or- 
der to prove that system (15)-(21) is exponentially stabilizable, let us first show 
that it is approximatively controllable. 

Theorem 5.1. System (15)-(21) is approximately controllable, in the space Hj!^ by 
boundary velocities ( G L^(0, oo; H'^/^(9<S)) n H-^(0, oo; H-^/^(5<S)) satisfying 

C • ndr = 0. 

ds 

Proof. Note that the operator formulation given by Proposition 4 does not enable 
us to write system (26)-(32) as an evolution equation. Instead of exploiting this 
operator formulation, we directly use the writing (26)-(32) and the definition of 
approximate controllability, as it is done in [17]. 

Let us show that if {uq, hi,ujo) = (0,0,0) then the reachable set R{T) at time 
T, when the control C describes L'^{0, oo; U^/'^{dS)) Ci R^{0, oo; U^/'^{dS)) with the 
compatibility condition 



dS 



is dense in the space ( J^) x 



C • ndP = 0, 

that we endow with the scalar product 



{{u, h' ,u!); {(j), k' ,r)) = / u ■ (j> + Mh' ■ k' + Iquj ■ r. 



J" 

For that, let (0^,fc'^,r^) be in R{T)^. We want to show that (0^,fc'^,r^) 
(0,0,0). 

Let us introduce the adjoint system 
d(f> 



dt 



— div (t((/), ip) = 0, 
div 6^0, 



in J" X (0,T), 
in J" X (0,T), 



(38) 
(39) 
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= 0, mdOx{0,T), (40) 

</) = fc'(t)+r(t) Ay, 2/e dS, t e {0,T), (41) 

- Mk"{t) = - I cr((/), V')ndr, t e (o,r), (42) 

JdS 

~hr'{t) = - [ yA a(0, i')ndr, t 6 (0, T), (43) 

JdS 

(t){y,T) = (f)^{y),yeT, k' (T) = k'^ (E M.^ , r(T) = r"^ e K^. (44) 
By integrations by parts from systems (26)-(32) and (38)-(44), we obtain 

u(T) ■ 0^ + Mh'{T) ■ fc'^ + Iou;(T) ■ r'^ = - f f (■ CT(<?i, V)ndrdt. 

J" Jo JdS 

Thus we deduce that if (cj)^ , k''^ , r'^ ) G R{T)^, then we have 

/ (•'^(^^"drdt = (45) 
Jas 

forallC6L2(0,oo;H3/2(a5))nHi(0,oo;Hi/2(95)) suchthat / C ■ "dr = 0. 

JdS 

This is equivalent to say that there exists a constant C{t) such that 

cr(0, ij)n = C{t)n on dS. (46) 

We can consider t/j = t/j — C instead of V', which does not modify the system (38)- 
(44). Thus we get cr((/), V')" = in L^{0,T]L^{dS)). Then equations (42) and (43) 
become 

k" = 0, r' = 0. (47) 

d(j) dijj 
Now, let us look at (pt = -rr- and -0* = with the condition 

ot at 

(T(0t,Vt)" = O on 95, (48) 

and the following homogeneous Dirichlet condition which is deduced from (40) and 
(47) 

0t = on aj" = 90 U dS. 

We can use an expansion of the solution of system 

-^-diva(0t,^O-O, in.Fx(0,r), 

div0t = O, inJ'x(0,r), 
0t=O, in9J"x(0,r), 

in terms of the eigenfunctions of the Stokes operator, similarly as it is done in [14]. 
Then the approximate controllability problem is reduced to showing that if 

—u/S.v + Vp = /iu in J^, 

div u = in 

V = {) on dF, 

(j{v,p)n = on dS, 
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with yu G M, then v = in Thus we get the following unique continuation result 
(sec [6]) 

4>t = in J^- 

Then we have only 

-div cr(0, V-) 0, in J" X (0, T), 

div(^ = 0, inJ'x(0,T), 

(^ = 0, inaC'x(0,T), 

= fc' + rAy, {y,t) e dS X {0,T). 

An energy estimate leads us to 

JJ^ JdS 

= k' ■ (7{(j),ip)ndT + r ■ / y A (T{<j>,tp)ndr. 
Jas JdS 

Combined to (46), we get 

J^\D{<P)fdy = 0, 

and thus -D(0) — 0. Using a result from [22, page 18], this implies 

(j) — k' + r A y in J-. 
The condition (40) enables us to conclude 

k' = 0, r = 0, and = in J". 
Then the proof is completed. □ 

Remark 2. The adjoint system introduced in (38)-(44) can be written in terms of 
operators; Indeed, by denoting <f) = (Fcj), k' ,r)'^ , we can formulate this system as 
follows 

-0' = A*<t>, 

ct>{T) = {n^,k'^yr, 

(Id-P)0 = (Id - P) (Lo(fc') + ^oW) , 

where A is self-adjoint, so that we can write 

ct>{t) = e(^-*)-^°(P0^,fc'^,r^)^. 

The adjoint operator of Z?o = ^Ba^o (whose expression is given in Proposition 4) can 
be expressed as ^?o0 = —a{z, 11)11^ where (2,7r) is defined as the solution of 

-i/Az + Vtt = (-Ao)0 in J", 
div z = in J^, 

z = ondF ^dO\JdS 

(see Lemma A. 4 of [16] for more details). Then, proving the result stated above by 
using the classical characterization of approximate controllability instead of using 
directly the definition would lead to other difficulties. 
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6. Stabilization and feedback operator. 

Theorem 6.1. For all X > and (uq, hi,LUo) G H;';^, there exists a boundary control 
C G L^(0,oo;H'^/-^((95))nH^(0,cxD;H^/^(35)) and a positive constant C - depending 
only on {uQ,hi,uJo) - such that the solution {u,h',uj) of system (26)-(32) satisfies 

/l', w)||h2,1(qo^)xH1(0,oo;R3)xH1(0,oo;R3) < C- 

Proof. Without loss of generality, we can choose A in the resolvent of A. Due 
to theorem 4.2, we know that the spectrum of —A is only a pointwise spectrum 
constituted of a countable number distinct eigenvalues, that we can order as follows 

3fiAi > 3?A2 > • • • > ?R:Xn > -A > SRA^v+i > • • • • 




Figure 1. The unstable modes of A are in finite number 
Moreover the generalized eigenspace of each eigenvalue is of finite dimension 
(see [9], Chapter III, Theorem 6.29, page 187). Let us denote by A(A,;) the real 
generalized eigenspace associated with A^ if A; £ R and with the pair (Ai,Ai) if 
3Ai 7^ 0, and let us set 

N oo 

H„ = 0A(A,), = A(A,). 

i=l i=N+l 

If E{\i) denotes the complex generalized eigenspace associated with A; and if 
{sj{\))i<j<m{\i) (where m(Ai) denoting the geometric multiplicity of A;) is a ba- 
sis of E{Xi), then A(Ai) is nothing else than the space generated by the family 
{Kej(Ai), Sej(Ai) | 1 < j < m{Xj)}. Let us observe that H„ is the unstable space 
of system (15)-(21) while is the stable space. Let us denote by P\ the projection 
onto the finite-dimensional unstable subspace H„ (parallel to the stable subspace 
Hs). If we project system (36) on H„, we obtain 

, / P^i \ / ¥u\ / Fu{;0) \ / uo \ 

-Px \ h' \ =APx\ h' + PxBoCo, Px h'iO) ] = Px \ h, . 

V - y V - / V-(0) / V-o/ 

(49) 

Due to Theorem 5.1, system (15)-(21) is approximately controllable in time T > 0. 
Thus the projected system (49) is also approximately controllable. Since it is of finite 
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dimension, it is also controllable. Let be a control such that Px{u, h' ,ijj){T) = 
(0, 0, 0), and let us still denote by its extension by zero to (T, cxd). Now, we notice 
that P\{Pu, h'jUj) is the solution of system (49) corresponding to Co if and only if 
P\{¥u, h' ,(2}) = e^*Px(P?i, /i', oj) is the solution of system 

^ / P^i \ /Fu\ fPui;0)\ fu„\ 

-Pa h' \=AxPx\ h' \+PxBxC, Px \ h'iO) \ = Px \ h, , 

U ; U ; U(o) J Uo/ 

(50) 

corresponding to the control C = e^*Co- 

Thus system (50) is stable. System (50) is the projection of system (15)-(21) onto 
its unstable space. Due to [23] and [13], system (15)"(21) is stabilizable by a control 
C G L2(0, oo; H3/2(a5)) n Hi(0, oo; U^/'^ldS)) if and only if its projection onto its 
finite-dimensional unstable subspace is stabilizable. The proof is complete. □ 

Remark 3. The regularity considered for the control C until now is L^(0, oo; H'^/^(3iS))n 
H^(0, oo; H-'^/^(9iS)). In the perspective of the nonlinear system in the second part 
of this work, we will actually need to consider more regularity, that is to say 

C G L2(0,oo;H5/2(c'5)) nHi(0,oo;Hi/2(a5)) 

So from now we will seek for the control ( in that space. 

Thus we know that system (50) is stabilizable by a control ( G L^(0,oo;S) 
satisfying 

/ C-^dr = 0, 

JdS 

where S is a subspace of finite dimension included in H^^^^dS) (see [2] for instance). 
So C lies also in H^/^((95). Moreover such a stabilizing control can be found in a 
feedback form. In other words, there exists 

ICx e £(H„;E!) 

such that the operator 

AxPx + BxJCxPx 

is stable on H„. Denoting U = (Pu, /i', oj), let us now consider the system 

PxV - {Ax + BxlCxPx)PxV- (51) 

When ICx is determined through an infinite time horizon control problem, by using 
the optimality system, it is easy to prove that the solution of (51) is actually such 
that /CaU € 11^(0, oo; S). Choosing ( = /Ca^aU in (36), this system becomes (51), 
and we have the estimate 

l|U|jH2,i((3(^)xHi(0,oo;R3)xHi(0,oo;R3) < C|| (^0, /^l , Wq) ||hi^ ■ 

We will still denote ICx the extension of ICxPx to the space H^, by setting 

ICxild-Px) = 0. 
The continuity of JCxPx provides 

||/CaPaU||hi(o.oo;H) < ^^||U||h„, 

and thus we also have 

II^AU||Hi(o,oo;HV2(a5)) < C[|U|| 
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Above all, such a choice changes the unstable operator Ax into A\ + B\IC\P\ which 
is stable. In the second part of this work, in fact we will only need the regularity 

C e L2(0,oo;H'5/2(5)) nHi(0,oo;Hi/2(5)). 

7. Definition of an admissible deformation from a boundary velocity. 

Given a control ( e L2(0, oo; U^/'^{S)) D R\0, oo; li'^/'^{S)) with the compatibility 
condition 

C • ndr = 0, 

dS 

(which can be chosen in a feedback form, as explained in the previous section), 
we want to define an internal deformation of the solid, denoted by X^, which is 
admissible in the sense of Definition 2.1. In particular we look for which satisfies 
the linearized constraints given by (22)-(24), and such that 
dX* 

e^'-gf{y,t) = ay,t), (y, t) e 55 x (0, oo). 

We also want the norm of e^*—g^ in the space 

L2(0, oo; n^iS) n Hi(0, oo; 'H.\S)) 
to be controlled by the norm of C in L'^{0,go;H^/'^{S)) nll^{0,oo;H^/'^{S)). 

7.1. Definition of a deformation satisfying the linearized constraints. We 

search for the deformation X^ by writing it as 



Xl{y,t) = y+ e-^^ip{y,s)ds, (52) 
Ja 

where the velocity (p{-,t) is the solution of the following elliptic system 

flip - 2div D{ip) = F{(p) in S, (53) 
ip = C on dS, (54) 

with 

F{p)iy,t) = ^ (^J^^2D{^)ndr^ + psl^' (^J^^y A 2D{^)ndT^ Ay. 

Proposition 6. For C e L'^{0,oo;U^/'^{S)) r\ll\0,oo;H^/'^{S)) satisfying 

C-nAT = 0, 



las 

system (53)-(54) admits a unique solution tp in L^(0, oo; H'^(5) H H^(0, oo; H"'^(iS)), 
for fj, > large enough. Moreover, there exists a positive constant C > such that 

ll<i5||L2(0,oo;H3(5)nHi(0,oo;Hi(5)) < C*!! CllL2(0,oo;H5/2(5))nHi (0,oo;Hi/2(5)) . (55) 

Besides, if Ci, (2 e L2(0, 00; H5/2(5)) nH^ (0, 00; Hi/2(5))^ and if pi and denote 
the solutions associated with Ci and C2 respectively, then 

\\^2 - </'l|lL2(0,oo;H3(5)nHi(0,oo;Hi(5)) < C*IIC2 - Cl llL2(0,oo;HV2(5))nHi (0,cx);Hi/2(5)) . 

(56) 

Proof. The proof of this proposition is given in section 7.2. □ 
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Remark 4. The compatibility condition assumed for the datum C is useless for the 
proof of Proposition 6, but contributes to making the mapping so obtained an 
admissible control (in the sense of Definition 2.1). 

Let us see that the mapping thus chosen is admissible. 

Corollary 1. For C e L2(0, oo; H5/2(5)) n Hi(0, oo; Hi/2(5)), the deformation X* 
provided by Proposition 6 and equation (52) is admissible for the linear system 
(15)-(21) in the sense of Definition 2.1. 

Proof. The constraints imposed in Definition 2.1 are equivalent to the following 
ones expressed in term of ip: 



Lp ■ ndr = 0, / ipdy = 0, y A tpdy = 0. 

dS Js Js 

Thus we have to verify that the mapping (p solution of (53)-(54) satisfies these 
constraints. The first constraint, which corresponds to (22), is satisfied thanks to 
the compatibility condition assumed for C. For the two other ones, let us recall that 
we have 



ydy = 0, 



5 



and since the tensor D{tp) is symmetric, we also have 

yAdivD{ip) = div {E>{y)D{ip)) . 
Then Equation (53) leads us to 

fi f ^dy = D{p)ndr) -2 f div{D{p))dy, 

Js \J dS / JS 

^I f yApdy = 2([ yhD{p>)ndT) -2 j diY {^{y)D{p)) dy, 

J S \J dS J J S 

and thus in using the divergence formula we get the two other constraints. □ 

7.2. Proof of Proposition 6. Instead of solving directly system (53)-(54), let us 
first consider a lifting of the nonhomogeneous Dirichlet condition. We set w the 
solution of the following Dirichlet problem 

div w = in iS, 

w = C on 95, 

with the classical estimates (see [7], the nonhomogeneous Dirichlet condition can 
be lifted by Theorem 3.4 of Chapter II, and the resolution made by using Exercise 
3.4 and Theorem 3.2 of Chapter III): 

I|w||h2(5) < C'IICIIh3/2(<s), 
I|w||h3(5) < C'IICIIhv2(5), 

Kt||Hi(5) < C'IICt||Hi/2(5). 

Then in setting ^ = — w, we are interested in solving the following system 
- 2div D{(\)) = F{<\)) - ^w + Aw + F(w) in 5, 

= on 95, 
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for some fj, > large enough, in the space H^(<S) in a first time. A solution of this 
system can be obtained as a fixed point of the following mapping 

N : H2(5) ^ H2(5) 

where </> is the solution of the classical elliptic system 

ficj) - 2div D{(p) = F{i') - fiw + 2div ^(w) + F(w) in S, (57) 

4:^0 on dS. (58) 

For proving that this mapping is well-defined, let us give some preliminary estimates. 
The equality 

2D{(t)) : D{(j)) ~V(t) -.Vet) = div ((0 • V)(/) - (div 0)</.) + (div c/))^ 
leads in Hj!)(iS) to 

l|V</.||^.(S) < 2\\DmlHS), 
and then the Poincare inequality provides a positive constant Cp such that 

We can estimate the norm H^(7^) as follows 

< Ci (||divi?(0)||2,(^^ + C^WDml^s)) ■ (59) 

The trace of D{ip) on dS which appears in the expression of F{ip) can be estimated 
as follows 

\\D{ip)n\\i^2i^QS) < C'2||V'||h3/2+-(5): 

with a = 1/2 + £, for some £ > which can be chosen small enough. Thus, by 
taking the inner product of the equality (57) by div £)((/)), we obtain 

^^\\Dmhis) + ndWD{<p)\\i,^s^ < c (ii^^(v)iil^(5) + iiw||h^(5)) wdivDm^^.^s) 

< c (p(V)^llL^(a5) + Mms)) + Iwd^Dmh 

By using (59) in the left-hand-side, it gives 

i2fi~3C'^)c^\\Dml.^s) + m\ims) < 2C1C3 {m^,^^s)\\Dml2^s) + Mims)) 

We now use the Young inequality on the right-hand-side, by introducing some 5 > 
which can be chosen as small as desired, as follows 

2CiC3p(v.)||^^(25")||v^||^V) ^ j\mms) + [^^) -WDmhis) 
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with p = 2/(1 + 2e) and q = 2/(1 - 2e). Then we have 

{2^,-3c'^)c^\\Dmhis) + m\H^is) < -iiv^iiH^(5) + (^-^j -WDmhis) 

+c||CIIh3/2(s)- 

Thus, by setting 

bh = {0 e h2(5) I (2a* - sc'p)c4Dmi.^s) + mh^(s) < r} > 

and by choosing S > small enough, and fi > and R > large enough, we can see 
that the ball B^j is stable by the mapping N. By the same inequalities we can see 
that N is a contraction in B^j, and thus N admits a unique fixed point in H^(5). 
The same method can be applied in order to prove the regularity in H'^(5). Indeed, 
since we have the equality V(divZ)(0)) = divZ?(V0), the gradient satisfies a similar 
equality, as follows 

l^i\7(p - 2divD{\7(p) = VF{(I>) - fiVw + 2divD{Vw) + VF(w) in S, 

so that we can show that V(f> lies in H^(iS). Then we have the estimate 

||0||h3(5) < C\\w\\h3(s), 

and since ip = (j) + w, we have 

||¥'||h3(5) < C'||w||h3(5) < C'IICIIhV2(5), 
ll'li'l|L2(0,oo:H3(5)) < C'||C||l2(0,oo:HV2 (5)) • 

The estimate which deals with the time-derivative of (j) can be obtained easily. 
Indeed, by taking the inner scalar product of the equality 

H(l)t - A(/)t = F{(f)t) + Wf 

by (jjt, we notice that the contribution of the right- hand-side force vanishes, as 
follows 



because (jjt satisfies the constraints 

0, / y A 0t = 0. 



By this means we get easily 

|I<<2||h1(0,oo;H1(5)) < C'IICIIh1(0,oo;H1/2(5)). 
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